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Abstract 

Charged particle energy spectra in e^e~ annihilation are compared 
with the analytical predictions from the QCD evolution equation in 
the Modified Leading Log Approximation. With the nonperturbative 
initial condition shifted down to threshold as suggested by the Local 
Parton Hadron Duality picture a good description of the data from 
the lowest up to highest available energies results. The two essential 
parameters in this approach are determined from a moment analysis. 
The sensitivity of the fit to the running of and to the number of ac- 
tive flavours (including a light gluino) is demonstrated. For very high 
energies the theory predicts a scaling behaviour in certain rescaled 
variables ( "^-scaling" ) . The data show an approximate behaviour of 
this type in the present energy range and come close to the predicted 
asymptotic scaling function for the small particle energies. 
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1 Introduction 



The perturbative QCD describes well the properties of hard processes. The 
problem of the soft limit of the perturbation theory and the transition to the 
hadronic final state, however, is still not solved at a fundamental level. The 
studies of the global properties of the hadronic final states gave support to 
the idea that the colour confinement mechanism should be rather soft . 
Meanwhile, the close similarity between hadronic and partonic final states 
has been found for a large variety of observables, provided the parton cascade 
in the perturbative calculation is evolved down to rather low virtualities of 
the order of a few hundred MeV; this phenomenon is called "Local Parton 
Hadron Duality" (LPHD) (for a recent review, see 

An interesting prediction concerns the momentum spectrum of the par- 
ticles in a jet at very high energies, namely the approximately Gaussian 
distribution in the variable ^ = ln(Pjet/phadron), the so-called "hump backed 
plateau" with a suppression of low momentum particles following 

from the soft gluon interference @J^. 

Within the LPHD picture, not only the very high energy behaviour of 
the spectrum is considered. Rather, one starts from the initial condition 
near threshold for the parton cascade where only one parton is present and 
derives the distribution at higher energies using the appropriate QCD evo- 
lution equation. In this case - apart from the overall normalization - the 
theory has only two essential parameters, the QCD scale A in the running 
coupling (at the one loop level) and a cut-off Qq in the transverse momentum 
of the gluon emission. In the Double Logarithmic Approximation (DLA) one 
takes into account the leading contributions from the coUinear and soft singu- 
larities which dominate at very high energies and determine the asymptotic 
behaviour of the observables. At present energies the next-to-leading correc- 
tions of relative order ^/al are important and are fully taken into account 
in the Modified Leading Logarithmic Approximation (MLLA). An explicit 
analytical expression for the particle energy distribution is found in this ap- 
proximation for the special case Qo = A, the so-called limiting spectrum 0. 
Remarkably, the experimental momentum spectra in the PETRA [^| and 
LEP |jTl| energy range gave support to the perturbative predictions 0, and 
this success has been continued recently to the higher energies at LEP-1.5 



g,0,|g and [|g) and the TEVATRON 1T§. 
The energy evolution of the spectrum and its connection to the initial 
condition at threshold is most conveniently studied using the moment repre- 
sentation of the spectrum. Recently, the first such moment analysis, based 
on the analytical results [jT^, has been performed over the full cms energy 
range available in e^e~ annihilation |T^. There are a number of advantages 
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of the moment analysis over the analysis of the spectrum itself: (a) there are 
analytical formulae which keep the two essential parameters of the theory 
Qq and A independent |T^, whereas explicit formulae for the spectrum are 
only available for the limiting case Qo = A; (b) the moments evolve with 
energy independently of each other. Their absolute size is determined by the 
initial condition at the threshold of the process; (c) the moments of order 
g > 1 are independent of the overall normalization and depend only on the 
two essential parameters of the LPHD approach; (d) the difference between 
the theory with running coupling and an artificial theory with fixed coupling 
can be studied directly and (e) the flavour thresholds can be included in the 
theoretical calculations. As the calculations of the moments include the soft 
part of the spectrum, the mass effects have to be taken into account. 

Having studied the behaviour of the spectrum towards low energies it is 
also interesting to study the limiting behaviour of the spectrum for very high 
energies as the theories typically predict simple asymptotic expressions or a 
characteristic scaling behaviour. Then one can investigate to what extent this 
behaviour appears already at available energies. The asymptotic behaviour 
of the QCD jet is obtained in the DLA. There is indeed a finite scaling limit 
in certain rescaled variables and the scaling function can be calculated . 

The purpose of this paper is twofold. First, we present more details of 
our previous moment analysis |T^ (see also |T^) as well as further results 



concerning all points (a) - (e) above. Secondly, we give the description of 
the shape of the spectrum and discuss the approach to the asymptotic limit 
in the appropriate scaling variables. 



2 Theoretical description of particle spectra 
2.1 Evolution equations 

The multiparticle properties of a QCD jet is conveniently described by the 
generating functional Qq, {u{k)}). The arguments are the hardness scale 
K of the jet, which is defined in terms of jet momentum P and opening angle 

eby 

K = 2Psm — ^ PQ (1) 

where the approximation holds for small angles; Qq is the lower cut-off in the 
transverse momentum of the emitted parton 

k± > Qo. (2) 



2 



The inclusive densities can be obtained by functional differentiation of Z over 
the probing functions refering to parton momentum k. Our analysis is 
based on the evolution equation for the gluon jet as discussed in |[T^. We 
rederive here the basic equations for the energy spectrum whereby we start 
from the evolution equation of the generating functional and show the various 
steps of the approximations. 

The evolution of the generating functional with the hardness scale k is 
given by the "Master Equation" [p|,pT| 

-IT— ^a{^,Qq) = \Y^ f dz 
am K 2 ^ JO 

$^(^) [Zb{zk, Qo) Zc{{l - z)k, Qo) - Za{k, Q,)] . (3) 
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This equation describes the evolution of the jet by the splitting of the primary 
parton A into partons B and C at reduced scales zti and (1 — z)k] ^^(-z) 
denote the DGLAP splitting functions. The integral over the secondary 
energy fraction z has to respect the limits (H) for the transverse momentum 
/cj_ = z{l — z)k. 

The Master Equation (3) yields results complete within the MLLA at 
high energies. Furthermore, the energy conservation constraints are taken 
into account. Because of (0) the evolution starts at k = Qo- Then the initial 
condition for solving this system of equations reads 

^A(/t,Qo; {m})U = Qo = UA{k = P) (4) 

which means that there is only one parton in the cascade at threshold. 

The solution of the evolution equation (3) together with the boundary 
condition yields perturbative expansions in in all orders, which for 
many observables can be resummed and exponentiate at high energies. The 
MLLA takes into full account the leading and next-to-leading contributions 
in the expansion of -Jal in the exponent at high energies. 

To this accuracy one can actually neglect the energy recoil and replace the 
1 — 2 by 1 in the arguments of Z and furthermore neglect the z-dependence 
of Z and of as{kj^{z)) in the regular z-integrals. With these simplifications 
one obtains a simplified form of (3), the MLLA master equation 

The evolution equation for the single inclusive density in the secondary 
parton energy k is obtained from Z by functional differentiation 

D{k,K,Qo) = j^^Z{k,Qo,{u})\u=i. (5) 
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After differentiation of the MLLA master equation |^ we obtain the 
evolution of the parton densities in quark and gluon jets 



dlnK 



Dg{k,K) = f —'yQ{zK)Dg{k,K) 

J X Z 



d 



dhiK 



Dq{k,K) = ^(^J^ ^-f^{zK)Dg{k,K) -^-f^{K)Dg{k,K)^ 



(6) 



Here x = k/P, Uf and Nc denote the number of flavours and colours respec- 

—1 1 

tively, Cp = ^'^'^ -^R ~ 2- '^^^ anomalous multiplicity dimension 70 is 

given in one-loop approximation by 

^°^^^) = — n — = bHkjAy ^ = y^^-r^- 

The initial conditions for these evolution equations follow from 

DA{k,K)\,=Q,=6{k-PA) , A = q,g (8) 
It is convenient to write these equations in logarithmic variables, namely 

>^ = ln7:-^ln— , ^ = ln- = ln— , A = ln— , 

Qo Qo X k A 

zP zPQ 
^' = In — = ^ + In = In — — = F + In ;z. (9) 

k Qo 

Then, after redefining D{k, P, 6) D{^, Y) and D{k, zP, 6) D{^', y') we 
obtain after integration of (|^) with (|DQ 

Dg{i, Y) = 6(0 + dC dy'^liy' + ODg{i', y' + O 

Dg{^, Y) = 5(0 + f f d^' dy'^liy' + ^ V + ^0 

-\[^ dy^l{y)Dg{i,y)\. 



""^We use here the same notation for both densities -D(fc) = dn/dk and D{£) — dn/d£^ 
kdn/dk. 
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The double integral terms in these equations originate from the singular 
parts of the splitting functions ~ 1/z and represent the leading double 
logarithmic terms of the DLA, the single integrals include the MLLA cor- 
rections from the finite parts of the splitting functions within the required 
approximation. In the DLA the spectra in quark and gluon jets are related 
by 

D,{^, Y) - 5(0 = ^(D.it Y) - m) (11) 

Replacing Dg by its leading order contribution j^Dg in the nonleading term 
of eq. (10) for Dg, we obtain the integral evolution equation for Dg 



Dg{^, Y) = 5(0 + d^' dy'^l{y' + ^D^'^ V + ^ 
Jo Jo 



Y 

dyil{y)Dg{i,y) 



(12) 



ANcJi 

with a = Y-^c + The corresponding differential equation reads 



^d^ ^ dYj 



dY 



For a = it corresponds to the DLA. In order to solve this equation one can 
introduce the Laplace transform {D = Dg) 

m,y)= ^e-^D^{Y), (14) 

where the integral runs parallel to the imaginary axis to the right of all 
singularities of the integrand in the complex a;-plane. Then D^{Y) fulfils the 
ordinary linear differential equation 

(^ + ^) ^DUY,X)-l'oiY)D^{Y,X) 

This is the basic equation used in fl^ to derive the parton distributions and 
their moments for gluon jets. For quark jets the high energy approximation 
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Dg = j^Dg is taken. The differential equation (13) has also been analyzed 



m 



One approach towards a solution of this equation is based on the anoma- 
lous dimension ansatz 



D^{Y, A) = D^{Yo, A) exp (^J^^ dy-f^[as{y + A)] 



(16) 



which yields a differential equation for with two solutions, one of which 
dominating at high energies. Alternatively, one can find the solution of (|TB|) 
directly in terms of hypergeometric functions [ll7|,3 . 



2.2 Moments of parton distributions 

There are various advantages of studying the moments of the parton distri- 
butions as discussed in the Introduction. Analytical predictions for QCD jets 
have been presented in [[7|,p!7|. 

2.2.1 Definitions and evolution equation 

The unnormalized moments A^^ of the ^-distribution are defined by 

Mg{Y,X)= r diim{i,YA) (17) 

JO 

with M.Q equal to the average parton multiplicity AT. They are closely related 
to the Laplace transform 

D^iY, A) = r die-^^ D{i, r. A), (18) 
Jo 

by 

Mg = {-ir^^D^{Y,X)l=,. (19) 
The moments < .^'^ > of the .^-distribution are then defined by 

<i\YA)>=^. (20) 

One also introduces the cumulant moments KqiY^X)] the moments of 
lowest order in q are given by 

=< e e, = a' =< (e - 0' >, 

^3=<(e-0'>, i^4=<(e-eT>-3a^ (21) 
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The reduced cumulants are defined by kg = Kg/ a'', in particular the skewness 
s = and the kurtosis k = k^^. The cumulants Kg of general order q can be 
derived from the expansion: 



\nD^{Y,X) = Y^Kg{Y,\) 

9=0 9 



[22] 



and therefore 



Kg{Y,\) 



d V 



(23) 



a;=0 



At high energies one term of the type (|T6|) dominates and one finds 
Kg{Y, A) = Kg{Q, A) + dy [-^ 7^[«.(l/ + A, rij)] 



(24) 



where the evolution starts at F = according to the LPHD picture. 

The moments can be obtained via eq. ([T9| ) from the evolution equation 
for with the appropriate boundary conditions at threshold. These are 
obtained from the Master Equation (3) with (^ and (|) as 



Ar(0) = 1, A^(0) = 
< ^9^0) > = < ^"'(0) > = for g > 1. 



(25) 



On the other hand, the boundary conditions implied by eqs. (10) and ([12|) 
are different in case of the multiplicity AT because of the approximations 
involved in this integral equation, whereas they remain the same for the 
higher moments. For A > we find 



Ar(o) 



Ar'(o) 



B 

T 



< 



(26) 



The corresponding result is also obtained for in [T^ . This yields a mini- 
mum of the multiplicity above threshold with an unphysical value Ar(y„) < 
1. In case of the limiting spectrum with A = both the position in Y and 
the value of the minimum of the multiplicity AT tends to zero and one ob- 
tains the boundary conditions ^"(0) = and AT (0) > Ofl^. These problems 
can be traced back to the approximations in the second integrals in (10) 
and (0): the correct limits (25) would be obtained, for example, using the 
modifications (e-terms) proposed in[JT . 



It should also be noted that the limit A ^ is only possible because 
of these approximations. Otherwise, whereas the energy dependence of the 
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multiplicity 77 ~ exp(cV^ + A) has a smooth limit for A ^ 0, the absolute 
normalization (the prefactor) would diverge for A ^ 0. This can be seen, for 
example, in the DLA from ( p2D with B = 0. It will be interesting to study 
the effect of these approximations in more detail. The subsequent analysis is 
based on the approximate eqs. ( |T2|) and (|l^) with the corresponding boundary 
conditions. 



2.2.2 Moments for running coupling 



The results obtained from the full solution of (|T^) can be written for arbitrary 
parameters Qo and A as: 

< e >= ^ E (^) {nA'' + N,Ri'^) (27) 

where A^i, iVa, L^'^ and Ri^^ are known functions of a, 6, F + A and A whose 
explicit expression depends on the order q [jl^ (see Appendix Al). For the 
special case of the limiting spectrum, where the two parameters Qo and 
A coincide (i.e. A = 0), the expressions simplify and all moments can be 



expressed in terms of the parameter B = a/b and the variable z = JlQNcY/b. 



The general result for the g-order moments is then the following ||T7| : 



^ = P^^\B + l,B + 2,z) + 2 + 1, 5 + 2, z) 



:28) 



where Pq"^^ Pi'^^ are polynomials of order 2(g — 1). Expanding the Bessel 
functions and Pi'^\ one obtains a series in 1/^/Y. The leading and next- 
to-leading order results in this expansion (see also ^) determine the high 
energy behaviour; the remaining part of the series however is still numerically 
sizeable at LEP energies (10% contribution to ^ and cr^) and increases towards 
lower energies. We therefore included the full result (^), also for p/'^'' (for the 
explicit expressions, see Appendix A2). The average multiplicity of partons 
is given in this approximation by: 

1~B 



JfLs = r(P) (^-j Ib+i{z). (29) 
2.2.3 Moments for fixed coupling as 

In this case the differential equations (13) and (0) can be solved exactly 



181. For D,,. one finds 



D{u;, Y) = [ ^'^^'^ smh{uY) + cosh(cDr) j e-("/2+^)^ (30) 
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where uj = y (a;/2 — r/)^ + 7q. Differentiation (|T9D yields the moments in the 
form 



fix 



cosh (70V) + — sinh (70!^) 
7o 



exp {-riY) (31) 



< 0» >= cosh (70^) + B, sinh (70^)] !^^Pp^ (32) 

A/ fix 



fix 

2 



7o = V7o+r ^ 



8iVc 47r 



The coefficients Ag, Bg are polynomials of order q in F, given for g < 4 in 
the Appendix A3. For 77 = 0, one obtains back the DLA results, especially 
J^fix = cosh(7oF). 

These results are obtained using the original boundary conditions, eq. (25). 
As pointed out above, from the MLLA equation (|1^) the boundary condi- 
tions ( p6l) are obtained. We have also studied the results following from 
these other boundary conditions: in this case the moments would be shifted 
towards larger Y at most by 1 unit and the overall description of the data 
would become worse. 

At high energies the results (32) greatly simplify and the moments read: 

< >^ ^^^^ (33) 

i.e., they behave like < >~ F^. For the cumulant moments the leading 
terms cancel and 

K,J^x ~ Y . (34) 

This can be easily seen by noting that at high energies eq. ( POD can be 
approximately written as 

D{uj, Y) ~ exp [7^y] with 7^ = -(| + r^) + (35) 
Since the coupling in this case is frozen, the anomalous dimension 7^ does 



not depend on y, the integral in eq. (p^ becomes trivial and eq. (|34D follows 
directly. 



3 Moment analysis 



9 



3.1 Determination of moments and mass effects 



In order to determine the moments one has to integrate the ^-spectra over 
the full range. Here one faces a problem for small momenta because of mass 
effects as already discussed in [113 . In the theoretical calculation the partons 



are treated as massless with k± > Qo, therefore ^ = ^£; = ln^<y = ln^, 
i.e. ^ has an upper limit. On the other hand, the experimental data usually 
refer to the distribution in particle momentum p or = In ^, and is not 
limited from above. 

For identified particles with known masses one can easily construct the 
energy distributions, but there is no reliable prediction yet for the mass 
dependence of the identified particle spectra within the theoretical framework 
considered here. Therefore we restrict ourselves to charged particles and we 
give them a common effective mass mh- If this mass is taken equal to the cut- 
off Qq, then for partons and charged particles there is the same upper limit 
C,E = y ■ The hadron spectra in the jet from the primary parton A are then 
calculated from the parton distribution according to the LPHD hypothesis 

from mm 



Eh — 3 = KhEp — = KhD\[e.E, Y) (36) 

dph dpp 



with Eh = \Jp'h + Qo and Kh an unknown normalization constant to be 
fitted by the data. Eq. (^) leads to the correct relation for p ~ ^ Qo, 
independent of the mass; for small momenta it yields a finite value for the 
invariant density E-^ in agreement with the data [|18| (see also the discussion 



on this point in ) . As Qo has the meaning of a transverse momentum cut- 
off for partons, it could be thought of for hadrons as an effective transverse 



mass TTiT = \J Tn^ + p\ , which is larger than the particle mass itself. 

The effect of introducing an effective particle mass on the shape of the 
spectrum is shown in Fig. (|ip. Here the inclusive charged particle .^-spectra 
with different mass assignments, namely Qo = 0, Qo = "^tt and Qo = 270 
MeV, are shown. The latter value is suggested from the moment analysis |jl8 



see below. The upper limits of ^ corresponding to the above effective masses 
Qo are also shown in the figures. As can be seen from these figures, the 
rescaling procedure is relevant in the soft region only, where the kinematical 
boundary becomes important. The effect is stronger at the lower cms energy 
of y/s = 14 GeV [0 as compared to ^/s = 91 GeV|]ll|. In the latter case 



the separation of curves in the measured range is hardly visible. 

The moments < ^'^ > are determined from the spectra Edn/dp vs. S,e 
after appropriate transformation of the measured Xp spectra and therefore 
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depend on the chosen effective mass mh = Qo- For the unmeasured interval 
near — Y (small momenta), we added an extra-point with coordinates 
{{y + ^iast)/'^, D{^iast)/'2±aD{^iast)} to linearly interpolate between the last 
measured point of coordinates {^last, D{C,iast) ± o'd^iast)} and the limit {Y, 
0} imposed by kinematics. The errors of the moments are determined from 
the errors of -D(^i) and the errors of the central values of in each bin, 
taken as half the bin-size. 

Correspondingly, we obtain the multiplicity AT^; as integral over of the 
full spectrum Edn/dp. As expected from Fig. (P, its difference to the usual 
particle multiplicity J\fch decreases with rising cms energy, from 30% at ^/s 
= 3 GeV to 10% at LEP energy. The MARK I data point at = 4.03 GeV 
shows an anomalous decrease of up to 50%, which may be partly related to 
charm thresholds effects. The result of this moment determination using the 
effective mass (5o=270 MeV is presented in Table 1; this value of Qo results 
from the fit discussed below. 



3.2 QCD description of moments for running as 

3.2.1 Determination of the parameters Qq and A 

The study of cumulant moments of the charged particle energy spectra al- 
lows for the first time the unconstrained determination of the two essential 
parameters which enter the theoretical predictions, namely Qo and A (or A) 
for the running as model and Qo and 70 for the fixed as model. Let us start 
with the running as case with results shown in Fig. (^, the analysis of the 
fixed as results, shown in Fig. (|]), runs parallel and will be discussed below. 

Fig. (0) shows the mean multiplicity, the average value C,e and the disper- 
sion cr^ extracted from the experimental data [|10|,pl],p!2tp!3|,p!4tp4|,| 



as a function of the cms energy for three different values of the parameter 
Qo] the theoretical predictions for the cumulant moments using the given Qo 
still depend on A (or A). The predictions are calculated for the number of 
flavours Uf = 3. For the particle multiplicity we use 

4 

J^E =Ci-2Afpart + C2 (37) 

y 

with arbitrary parameters Cj and parton multiplicity J^fpart from eq. (|29|) (the 
factor 4/9 is for the quark jet and 2 is for the two hemispheres). The two 
normalization parameters are determined to let the curve go through the 
lowest and the highest energy data points. The parameter ci corresponds 
to the Kh factor in eq. (|36|) , whereas the additional parameter C2 has been 
introduced to allow for a finite multiplicity at threshold as in (25). It is 
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important to note that the higher moments (g > 1) describe the shape and 
do not depend on the normahzation, so they are unaffected both by the 
systematic experimental uncertainties of the overall normalization and by 
the theoretical uncertainties associated with the factor. 

The mean multiplicity data in Fig. for each chosen Qq, can be prop- 
erly described by the theoretical predictions for any value of A in the range 
< A < 4. Let us stress that the theoretical predictions at parton level 
strongly depend on A, but one can obtain in all cases good fits of the exper- 
imental data by adjusting the free parameters Cj in eq. (^). Therefore the 
mean multiplicity data alone do not determine the parameter A. Looking 
now at the first moment ^, we observe that for each chosen Qq a suitable 
value of A can be found which provides a good description of the data. In 
order to fix both parameters, one has then to include the dispersion cr^. As 
can be seen from the figure, a lowering of Qq shifts both the and the cr^ 
data downwards. On the other hand, an increase of A yields lower values for 

but larger values for a^, as one can see by inspection of the first two terms 
in the expansion of and cr^ in A [IS]. 



The parameters Qo and A are determined from a x^-minimization. To 
compute the we have used the first 4 moments (g > 1) of the inclusive 
energy spectra for charged particles, but not the mean multiplicity because 
of its larger systematic errors and the need of two more normalization pa- 
rameters for its theoretical description. The minimum of is obtained for 
the limiting spectrum (A ^ 0), and the parameters are estimated as: 

Qo ^ A ^ 270 MeV. (38) 

The minimum value of the /d.o.f., neglecting the correlations among 
the cumulants of different order, is found to be 1.8 (with about 70 d.o.f.). So 
we do not obtain a "perfect" fit of the data, but considering the small errors 
of the lowest order moments and the small number of parameters, a very 
satisfactory description of all moments over a large energy range is obtained. 
In view of the systematic uncertainties of the fit we estimate the errors of 
the parameters from the limiting case in which the theoretical curves miss all 
data points by about one standard deviation. This yields the conservative 
estimate 

AQo ^ AA ^ 20 MeV. (39) 

Alternatively, this result can also be transformed into a limit on A = ln((5o/A) 

A < 0.1. (40) 

Our result (|38|) is slightly larger than the ~ 250 MeV obtained by the OPAL 
Collaboration |TT[ from a fit of the measured ^ distribution to the limiting 
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spectrum adjusting only one parameter. The small difference results from 
the inclusion of lower energy data in our fit. 



3.2.2 Discussion of the fits 

The cumulant moments up to order q = 4 together with the correspond- 
ing predictions for the limiting spectrum with Qq = 270 MeV and nj = 3 
are shown in Fig. (0) by the solid lines. The values of the moments and 
the corresponding theoretical predictions are given also in Table 1. A very 
satisfactory description is obtained in the full cms energy range available; 
small deviations in the first moment at very low energies are visible. This 
may signal some limitations of the approximations involved, in particular, 
the simplified relation between quark- and gluon-jets, in this region. 

It should be noted that the moments with g > 1 are determined by 
two parameters only which actually almost coincide. This should be con- 
trasted with the more conventional applications of perturbative QCD to the 
particle or parton spectra. There one starts at a finite energy Yq with a non- 
perturbative input distribution which in general requires a set of unknown 
parameters and then evolves this distribution to higher energies according to 
the predictions of perturbative QCD. In terms of moments this would require 
one adjustable value Kq{Yo) for each moment at the initial energy lo- 
in the application of LPHD one assumes the validity of the perturbative 
formulae for the moments down to small energy scales of order few hundred 
MeV, actually down to the threshold energy Qq, where the distribution func- 
tion is known to be simply the 5-function (10). In this limit all the higher 
moments are determined to be zero. Therefore the compact description with 
only two parameters is a direct consequence of the assumption that the the- 
oretical description can be continued down to these low energies. A free 
adjustment of the moments at higher energy Yq with a vertical shift of the 
curves in Fig. (^) would not improve the predictions essentially; the limiting 
spectrum with the absolute normalization at threshold gives indeed the best 
results. It should be noted that the previous applications of the limiting 
spectrum to fit the ^-spectra (for example |§|,|n|) rely on the same assump- 



tion that the QCD evolution can be continued down to the low scales Qo of 
few hundreds MeV where the initial condition is introduced. 

3.2.3 Flavour dependence 

The theoretical predictions shown so far were obtained with 3 active flavours. 
A possible source of uncertainty in the theoretical formulae is the number of 
active flavours to be used (see also |jIB|,p!5|). In the predictions for cumulants 
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(|2^ ) the number of flavours enters essentially through the running coupling 
asiy, rif). We neglect in the present discussion the additional explicit depen- 
dence at the percent level, which comes in through the parameter a at the 
next-to-leading order of the MLLA. 

The moments evolve at low energy with 3 active flavours and with 4 and 5 
flavours after passing the respective thresholds. The simplest approach would 
be to put the thresholds at the heavy quark masses, i.e., to increase n/ by 
one at ^ = niQ where y/s = 2P. However, let us recall that the argument 
of tts is the transverse momentum k± and kinematics forces k± < This 

suggests moving the thresholds to ^ = 4mQ (or towards even larger values, 
see for example pO|). 

In Fig. (^) we show the predictions from the limiting spectrum with the in- 
clusion of heavy flavours at the corresponding thresholds ^ = Amg (dashed 
lines). Above a heavy quark threshold the moments evolve according to (p^) 
with the respective number of active flavours Uf and match continuously to 
the moments with Uf — 1 active flavours below the threshold, in complete 




analogy to the inverse coupling constant l/as|^T|. We then write for the 
cumulants 

Here K^\Pjet) refers to the moment calculated with i flavours from thresh- 
old Qo up to jet momentum Pjet- The comparison of the upper two curves 
in Fig. (^ shows that the inclusion of the heavy quark thresholds does not 
modify dramatically the behaviour of the moments. Therefore, a reason- 
able approximation of the experimental data in the present energy range is 
obtained by taking into account only three active flavours throughout the 
full energy range. This behaviour of the theoretical predictions can be eas- 
ily understood from the representation of the moments (0) in terms of the 
anomalous dimension. Just above threshold the contribution of the new 
flavour to the y-integral is negligible. Most primary gluons are still emitted 
at a smaller scale below the new flavour threshold. On the other hand, an 
approximation with flve flavours in the full energy range would be a complete 
failure. 

We have neglected here the differences in the light and heavy quark frag- 
mentation. Whereas such effects occur in the fragmentation region, they 
are expected to be small in our application, as the soft gluon radiation is 
universal. 
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3.2.4 The effect of a light gluino 



By exploiting the flavour dependence one can also extract new information 
on the possible presence of additional light particles. There has been con- 
siderable interest in the last years in the question whether there is a super- 
symmetric gluino with a small mass. In a recent summary the gluinos in 
the mass range l|-3| GeV are considered as absolutely excluded, whereas 
lighter gluinos are allowed, except for certain ranges of lifetime. In a recent 
study of jet rates and jet angular distributions in the reaction e~^e~ — 4 
jets, such a possibility has been severely restricted, however. 

A sensitive probe of the presence of light gluinos is the running of as ||34|| , 
as each gluino changes the number of active flavours by 3. Here we show the 
sensitivity of the moments to the presence of a light gluino with mass around 
1 GeV. It is assumed that the effect of the light gluino comes in only through 
the running coupling and not through its effect on the final state structure. 
This can be justified by noting that gluino pair production - like quark pair 
production - does not contribute to the cascade evolution in leading double 
log order. In Fig. the lowest curve represents the predictions for the 
moments assuming the presence of one light gluino with 1 GeV mass, i.e., 3 
additional flavours at ^ > 4 GeV. 

The multiplicity AT can be fitted again by readjusting the normalization 
parameters in agreement with previous findings |^ . On the other hand, the 
predictions for the higher moments, especially with q = 2 and g = 3, are 
far off the data. The energy dependence of the moments in presence of the 
light gluino is weaker in the same way as the running of is weaker, as it is 
expected from eq. p^). We conclude that the existence of a light gluino is not 
supported by our analysis. However, it seems premature to definitely exclude 
such a particle at present from this study. There are some simplifications in 
the present analysis and our QCD fit without gluino is not perfect in a 
sense. Since the moments are very sensitive to the existence of a light gluino, 
a meaningful statistical test could be performed after a further improvement 
of the theoretical description. 



3.3 Results for fixed 



a. 



To see the effect of the running coupling in the inclusive energy spectra, let 
us now consider for comparison the corresponding model with fixed coupling. 
Fig. contains the same data as Fig. (^, but the theoretical predictions 
refer now to the MLLA with fixed as ll^]- The curves correspond to three 
different values of the coupling a^, i.e., of the anomalous dimension 70. For 
the multiplicity we take eq. (^) with TJpart given by eq. ([31|). Again, the 
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parameter 70 cannot be extracted from the study of the mean multiphcity 
alone. Including the other two cumulants, one can reproduce at best, choos- 
ing 7o = 0.64 (i.e., = 0.21), the multiphcity data and the energy slope of 
the first moment ^e, though not its absolute value. 

In Fig. d^) we show the predictions of the fixed-a^ model with 70 = 
0.64 (dashed lines), where the absolute normalization is determined again at 
threshold as in the case of the running coupling, see (32). We also investi- 
gated whether the agreement with data can be improved if the normalization 
at threshold is abandoned and shifted to a higher energy. To this end we 
introduced an additional parameter for each cumulant, which allows for ver- 
tical shifts of the curves; they have been chosen to fit the experimental points 
at ^/s = 44 GeV. In the following we will refer to this model as the shifted 
fixed-tts model. 

Contrary to the case of running as, the vertical shifts can improve the 
description of the moments for the model with fixed coupling, but only in a 
limited range of cms energies. Especially, the moments with q > 3 show a 
rather different trend with energy in comparison with the data. 

In an alternative investigation of the relevance of the running coupling for 
the inclusive energy spectra, it has been proposed to look at the predictions 
of Monte Carlo programs with and without running coupling |36|. It was 
found that the JETSET Monte Carlo with the standard hadronization 
phase but the coupling frozen at the value of = 0.2 in the perturbative 
phase describes the experimental data reasonably well throughout the PE- 
TRA/PEP energy range and deviations occur only at higher energies. At 
first sight this result seems to contradict our findings. Note, however, that 
in the JETSET Monte Carlo the perturbative evolution stops at a cut-off 
value of about 1 GeV, when the string fragmentation takes over. In our 
perturbative approach, we allow on the contrary the perturbative cascade to 
evolve down to the much smaller cut-off Qo ~ 270 MeV. It is in this low 
energy domain that the variation of the coupling is most pronounced. The 
running coupling becomes large especially for small k± so that particles tend 
to be produced coUimated. The perturbative calculations at low scales with 
running coupling seem to simulate the production and decay of resonances 
implemented in Monte Carlo programs like JETSET. This supports the idea 
that the parton cascade with running coupling down to small scales is dual to 
the cascade with a shorter perturbative phase but with hadronic resonances 



in the last stage 18 
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3.4 Rescaled cumulants 



In addition to the standard moment analysis performed in the previous sub- 
section, let us also consider the rescaled cumulants Kq/^. These quantities 
become energy independent in case of fixed coupling at high energies as fol- 
lows directly from (0). In particular for the first three rescaled cumulants, 
one has 

i 27^r/ + 7o ' i 47^r7 + 7o ' i 87^ 77 + 70 ^^^-^ 

Therefore these ratios exhibit more directly the differences to the case of 
running coupling. In Fig. the experimental data on these ratios, as derived 
from our moment results in Table 1, are compared to the MLLA predictions 
with running a^, then with the fixed and the shifted fixed-Og models. Once 
again, a good description of data is given by the MLLA model with running 
coupling. The fixed-a^ model shows the expected behaviour, i.e., the rescaled 
cumulants tend to a constant value at large cms energies; its predictions lie 
far away from the experimental data. The odd moments, which vanish in the 
DLA, approach their asymptotic limits more slowly. In the shifted fixed-a^ 
model previously described, the behaviour of the rescaled cumulants changes 
and the predictions become closer to the data. However, the deviations from 
the running predictions become obvious for the smaller and also the very 
high energies with F > 6, where the ratios reflect the different asymptotic 
trends. For example, the ratio K^/^ changes curvature when going from fixed 
to running a^. Data at these higher energies are becoming now available at 



the TEVATR0N||16|| and could give new information. 



4 Analysis of the shape 

4.1 Energy evolution of the shape 

The moment analysis has selected the solution with similar values for the 
parameters Qq and A; let us now consider the predictions for the shape of 
the spectrum itself. Fig. (|^ shows the inclusive energy spectra Edn/dp as a 



function of $,e, extracted from experimental data using the 

fitted cut-off parameter Qq = 270 MeV as the effective mass in the calcula- 
tion of the particle energy. The curves show the predictions of the limiting 
spectrum with the same value of the Qq parameter. The normalization has 
been fixed by choosing the integral of the spectrum to be equal to the average 
multiplicity according to the formula (^), the respective numbers are also 
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given in Table 1. The fit describes well the main features of the data in the 
wide range of cms energies 7 < ^/s < 140 GeV, especially in the region with 
smaller than the peak position. 

Some deviations of the fitted curves from the data can be seen for larger 
^^;'s, i.e., for smaller particle energies. At low cms energies the curves fall 
somewhat above the data near the peak position and below the data near 
the kinematic limit ^ Y)- This behaviour may be related to the fact that 
the limiting spectrum approaches a constant value and not zero for ^ — »• F, 
as expected for the exact solution of the evolution equations (10) or (IT^) . 

A iterative approximate solution of the MLLA equations which is valid 
in the soft region and goes to zero for —>■ Y has been given in |23[|: 



Y, X)mlla = Di^, F, A)|^L^exp 



Y 0/2 



7|M 



dy 



c 



(43) 



where 



D{tYA)\DLA = ^\n{l + ^^ 



1 + 



b 



ln(l + ^) 



In Fig. (H) the theoretical predictions from this approximation are compared 
with the same data as in Fig. (|^) for low particle energies E < 1 GeV. A 
rather good description of data in the soft region is obtained in this way, in 
agreement with previous findings on the quantity Edn/d^p ||23| . 



4.2 High energy limit of the spectrum and (^-scahng 

In QCD the asymptotic behaviour of the soft production phenomena can 
be derived within the DLA, which takes into account the leading collinear 
and soft singularities originating from the Bremsstrahlung processes. A well 

of the particle 



known scaling law of this type is the KNO scaling [^^ ,^ 
multiplicity distribution for rescaled probabilities and multiplicities. 

Within QCD the momentum spectra do not scale asymptotically in the 
Bjorken or Feynman variables x = p/P because of the Bremsstrahlung emis- 
sions with large transverse momenta, as is well known. Instead, they ap- 
proach a finite scaling limit in certain rescaled logarithmic variables. This 
proposal, put forward already more than ten years ago [^], has never been 
studied since. A scaling behaviour of similar type for angular correlations 
has been recently proposed |^ and found some support by experimental 
data El. 
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The asymptotic behaviour of the inclusive energy spectrum can be derived 
from the DLA evolution equation, i.e. ([T^) with a = 0. One finds a scaling 
law in the rescaled logarithmic variable ( for the rescaled spectrum 



^ Y 



(45) 



independent of Y. The logarithmic variables naturally occur in the asymp- 
totic limit as they absorb the Bremsstrahlung singularities in the evolution 
equation (P). 

The scaling function F{Q takes in DLA the following form pQ|: 



sinh /i 



with 



2C-1 



sinh 2/i — 2/i 
2 sinh^ /i 



(46) 



The function F{() is symmetric around ( 
of 1, and it goes as 



1/2, where it has a maximum 

Clnl/C for C ^ 0. In the MLLA model with fixed 
coupling explicit analytical expressions for the spectrum and for the average 
multiplicity are available [|l^ and the scaling function is given by F{() = 

s^ai-c). 

For further illustration of the scaling law (^), we recall the Gaussian 
approximation of the DLA [|,|| 



my) 



((r + A)3/2-A3/2)l 



exp 



(Y + A)3/2 - A3/2 



(47) 



As the multiplicity J^(Y) ~ exp( 
finds in exponential accuracy 



16Nc{Y + X)/b) at high energies, one 



2^ 2^ 



(4J 



which depends only on ( and not on Y. In case of fixed as the Gaussian is 
slightly narrower F{Q ~ 1 — 2(C — |)^. The DLA asymptotic scaling functions 
both for running and fixed ag and their Gaussian approximations are shown 
in Fig. (P). We notice that the two asymptotic curves for running and fixed 
Os do not differ very much quantitatively in the region around the maximum 
( ~ 1/2, but they do in the soft region, where the effect of the running of 
as becomes relevant p3[ and the corresponding curve shows indeed a steeper 



slope than the fixed-Os one. Notice also that the Gaussian approximation 
considerably deteriorates towards the limits C ~ and C ~ 1- 



It should be noted that only the transformed observable ( ^5]) and not the 
spectrum D{(, Y) itself approaches a finite limit. The spectrum D{(, Y) can 
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be approximated by a Gaussian in C with width cr^ ~ and therefore 

one finds asymptotically for D and its moments 

D{CY)^5{C-\) , <C>^Yi ^^^^ 



The MLLA results for the moments (p^)-(|66|) indeed approach this limit. 

In Fig. ( p!OD we show the experimental data for the observable (|^) at 
three different cms energies In this figure ^ is chosen to be and 

Y = \n{^/s/2Qo) with Qq = 270 MeV. Experimental data for the average 
charged multiplicity have been taken from p5| , ^|l4| . Fig. (|ll]) shows the 
same observable, but this time as a function of = C,e, where the energy is 
calculated with the effective mass Qo = 270 MeV. In this way, experimental 
spectra have a common kinematical boundary at C^; = 1 and can be compared 
with the theoretical predictions as in Fig. (J^). Fig. (p!T| ) also shows the 
theoretical predictions of the Limiting Spectrum normalized by the average 
multiplicity p7| ) and (^) at the same cms energies of the experimental data, 
as well as the asymptotic DLA prediction of eq. (|46|). 

One can see from these figures that the original energy evolution of the 
spectrum visible in Fig. (0) is largely removed if the scaling variables (|45|) 
are used. Some scaling violation remains, especially in the small ( region. 
The limiting spectrum reproduces well the small scaling violations shown 
by the data. It approaches the asymptotic limit very slowly and only at 
unphysically large energie^. The data at the present cms energies lie far 
away from the asymptotic DLA curve in the small ( domain; however, it is 
interesting to notice that the data for soft particles (p ~ 0, ^ ~ F, C ~ 1) 
are already close to the asymptotic limit. This result lends further support 
to the idea that very soft particles are basically determined by the nearby 
Bremsstrahlung singularities and not affected by the non singular terms in 
the splitting functions taken into account in MLLA nor by the recoil effects 
for a previous discussion of this point). 



see 



The energy dependence of the position of the maximum of the ( distri- 
bution, (*, is shown in Fig. (0) up to LEP-1.5 energies. Also in this case, 
the data closely follow the prediction of the limiting spectrum, which very 
slowly approaches the asymptotic DLA limit = |. 

One concludes that the data show an approximate scaling law in the 
presently available energy range. However, this scaling is preasymptotic, 
i.e., the asymptotic shape of the distribution is quite different from the one 
observed at present energies. Similar results on the existence of a large 



^For instance the scaling function F{() at C = 1/2, whose asymptotic value is 1, reaches 
0.8 at r - 68 (^/s - 1028 GeV) and 0.9 at F - 460 - 10^°° GeV). 
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preasymptotic scaling regime have also been predicted for the multiplicity 
scaling (KNO) and its violation [^,|^; in this case significant deviations 
from scaling are expected only around ^/s ~ 1 TeV ^ . 



5 Conclusions 

The perturbative QCD approach has been shown to describe well experi- 
mental data on charged particle inclusive energy spectra in e~^e~ annihila- 
tion. The new features of the moment analysis of the spectrum have been 
discussed. 

The first determination of the two independent essential parameters of 
the theory has been performed, whereby the best description of the data is 
obtained for Qq — A ^ 270 MeV with an uncertainty of about 20 MeV. The 
dependence of the moments on the initial conditions has been studied; if one 
takes the nonperturbative initial condition of the perturbative QCD evolution 
at the threshold of the process, a good description of the moments also in 
their absolute normalization is obtained in the full energy range available 
including the low energies of a few GeV. The sensitivity of the moments to the 
running of the coupling has been established by comparing the predictions of 
the full model with the predictions of a model with frozen coupling. The 
effect of heavy quark thresholds in the running of the coupling has also 
been discussed; a good phenomenological description of the data at present 
energies is obtained already by including only three active fiavours. The 
moments are very sensitive to the presence of a light gluino with mass around 
1 GeV, but there is no evidence for the type of effect expected. 

The inclusive energy spectra themselves have also been studied. The 
limiting spectrum with Qq = A is found to provide a good overall description 
of the data in a large cms energy range. An improvement in the soft region 
can be obtained by applying an approximation which takes into account the 
boundary conditions at ^ = F explicitly. 

Data in the presently available cms energy range support an approximate 
scaling law, the ^-scaling, predicted at asymptotic energies. Some violations 
of the scaling behaviour in this preasymptotic energy range are observed, as 
expected from the MLLA. The data in the soft energy region {( < 1) at 
available energies are already close to the asymptotic predictions. 

There are some simplifications both in the phenomenological analysis 
(for example, neglect of difference in light and heavy quark fragmentation) 
and the theoretical description. The respective improvements would yield 
a more stringent test of the theoretical scheme and possibly a quantitative 
improvement of the fit to the data. As the main result of the present analysis 
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we consider the good description of the data down to low scales - even if not 
entirely quantitatively everywhere - in agreement with QCD with running 
as and the LPHD picture. 
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A Explicit Formulae for Moments 
A.l Moments in the general case 

The functions entering eq. ( P?] ) are defined as follows[|T^: J\f = Afi + A/2 is 
composed of two terms increasing and decreasing with energy, respectively, 

Afi = T{B)T{l-B)z, Ib+i{zi)I-b{z2) (50) 

M = T{-B)T{B + l)zi (j^"" I-b-i{zi)Ib{z2) (51) 



z, = Vl6iV,(F + X)/b; Z2 = Vl6iVeA/6. 
The two terms build up the total multiplicity 

AT = (^^y [Ib+i{zi)Kb{z2) + Kb+i{z)Ib{z2)] (52) 

ov.^ 0(9) 



The functions and in ( [27| ) are computed from 

Lk^ = D,-k{B + l,B + 2, z^)Dk{-B, 1 - B, Z2) (53) 
Rk^ = ^.-fc(0, -B, ^i)Z^fe(0, B + l, Z2) (54) 

Z^.(^, c, z) = P(^)(^7, c, z) + -/^A^^)(^7, c, z) 
where P/'^'' are the polynomials 



Po^^^(^, c, .) = E «S1 ( 7 ) , P^\9. c,z) = Y. 



fc=0 fc=0 



(55) 



The coefficients of highest order are given by 



= pi^) = ±(B + ^) (56) 

whereas in lower order they can be found by solving a q x q linear system 
of equations. For the first four moments one finds explicitly ((a)n = T{a + 
n)/T{a) = a{a + 1) . . . (a + n — 1)): 



1-c' (c-2)2 (C-1)2 

2(c-3)3^(c-2)3 2(c-l)3 ^ ' 
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in terms of the expressions 

$(1) 



= -{n - 1)2 + gn 



(59) 



$(2) = 1 _ 2)3 + + - 2)3 + {g)2{n - 1)2 



(60) 



+ {lg' + lg + l)in-3), + igUn-2), 



(61) 



^^'^ = hn-7)s + ilg + l)in-6h + i^g' + ^g + f){n-5), 
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24. 



+ {2g' + 10(7^ + 14(7 + - 4)5 + (n - 8)4. 



(62) 



A. 2 Limiting spectrum 

In the special case Qq = A we find fl^ for the moments q < 4: exphcitly: 

^ 1 , BIb+2{z) 



Y 2 z Ib+i(z) 



(63) 



<e^> ^1 BjB + l) (3 + I) / _ 2B{B + 2) \ Ib^z) 



Y2 



Ib+i(z) 



(64) 



<^^> _ 1^35(5 + 1) 2S2(5 + l)(5 + 3) 



+ 



2^2 ^4 

3^ + 2 _ B{B + 1){B + 3) ^ 4B{B)i 



(65) 



^5+1(2;) 
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and 



<e> 

y4 



+ 



16 

2 

z 



+ a 



(4) 



2x 2 



(4) 



6 

1 + 



zj 



2\2 



where 



(4) 



= 10 



(4) 



l^(^ + l)(E + 3) 
4 1 



and 



(S)4 ' ^ (5 + 1)2 



^(4)^_6^+l.(i^)4 



(5-1)5 45-1 



" (5)^ ' - (5TI) 

and 

= —B^ + -B^ + 3.458335^ + 7.75^ 
+ 8.395835^ + + —B^ - -B. 

A. 3 Fixed coupling 

The first coefficients Ag, Bq in (32) read (here p = rj/jo) 

A,^{p' + 1)^ , Si = (-p2-l)^ + py ; 
A2 = (-3p^-p)^ + (3p^ + l)^, 

= (3p^ + p)^ + (-3p^ - + + 
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(72) 



A4 = (-105p^ + 30p3 + 27p)-^ + (75p^ - ISp^ -9)^ 

I670 I670 

+ 2(-5p^ - Up' + 3p)^ + {5p' + lOp' + 1)^ , 

1 Y 
= (105p^ - 30p3 - 27p):^^ + (-75p^ + 18p^ + 9) 



I670' 



167(1 



y2 y3 

+ (45p^ + 18p3 - 15p)— 2 + 2{-15p' - 2p2 + 1) 



I670' 



I670 
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Table Caption 



Tab. 1: The average multiplicity 71 e, the average value ^e, the dispersion 
(T^, the skewness s and the kurtosis k of charged particle inclusive energy 
spectra Edn/dp vs. for Qq = 270 MeV at various cms energies y/s. 
In brackets theoretical predictions of the limiting spectrum of MLLA with 
running a^; the second entry in the average multiplicity column contains 
the results of eq. (^) with eq. (^9]) ; the first one the results with C2 = 0. 
Errors on the average multiplicity data points include both statistical and 
systematic errors. Results at LEP and LEP-1.5 cms energies from 
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Figure Captions 



Fig. 1: a) Comparison of charged particle inclusive single particle spectra, 
Edn/dp vs. ^, at ^/s = 14 GeVfl^ for different mass assigment: the inclusive 
momentum spectrum pdn/dp vs. (diamonds) and the rescaled spectra 
Edn/dp vs. C,e, with i?^ = + Ql, Qq = 138 MeV (triangles) and Qq = 
270 MeV (squares). Also shown are the upper limits of C,e given hj Y = 
ln(-ys/2Qo); b): same as in a), but at ^/s = 91 GeV [p]] . 

Fig. 2: Dependence on Y of the three lowest order moments of the inclusive 
energy spectrum. The data of the mean multiplicity 77 e, average value C,e 
and dispersion of inclusive energy spectra for Qq = 138 MeV (triangles), 
Qo = 270 MeV (diamonds) and Qq = 350 MeV (squares) are compared with 
theoretical predictions of MLLA with running coupling with A = (limiting 
spectrum) (sohd line), A = 0.5 (dashed line), A = 4 (dotted line). The 
theoretical predictions for the mean multiplicity are computed from eqs. (|37|) 
and (H). 

Fig. 3: Same data as in Fig. 2 for different Qo parameters; comparison with 
theoretical predictions of MLLA with fixed coupling with 70 = 0.64 (solid 
line), 7o = 0.4 (dashed line), 70 = 1 (dotted line). The theoretical predictions 
for the mean multiplicity are computed from eqs. (^) and (pi]). 

Fig. 4: The average multiplicity 77 e and the first four cumulants of charged 
particle energy spectra Edn/dp vs. ^e, are shown as a function of Y = 
ln(A/i/2Qo) for Qo = 270 MeV at various cms energies (see Table 1). The 
curves show the predictions of the limiting spectrum with Qq = 270 MeV 
with 3 active flavours (solid line), the number of flavours Uf variable with 
the heavy quark thresholds at y/s/2 = Amq (dashed line) and the inclusion 
of a light gluino with mass of 1 GeV (dotted line). In all cases, the mean 
multiplicity is computed from eqs. ( pT]) and (pOl) . 

Fig. 5: Same data as in Fig. 4, but the curves show the predictions of the 
limiting spectrum (i.e. Qo = A) of MLLA with running as (solid line), of 
MLLA with fixed (dashed line) and of the shifted fixed model (dotted 
line); nj = 3 everywhere. 

Fig. 6: Rescaled cumulants Kg/^^; as a function of y = ln(A/i/2(5o) for Qo = 
270 MeV. Data as in Fig. 4 are compared with the corresponding predictions 
of the limiting spectrum of MLLA with running as (solid line). The dashed 
curve shows the predictions of the model with fixed as (70 = 0.64); in this 
case, the rescaled cumulants approach constant values at high energies (0.534, 
0.344, 0.528 respectively, with the chosen value of 70). Predictions of the 
shifted fixed as model are also shown (dotted line). 
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Fig. 7: Comparison of inclusive spectra Edn/dp vs. ^e, with E"^ = p"^ + 
Qo^ Qo = 270 MeV at various cms energies with predictions of the limiting 
spectrum {Kh fixed at each energy from the fit of eqs. ( |57D and (^)). Each 
curve is shifted up by 0.5 for clarity. 

Fig. 8: Comparison of inclusive spectra Edn/dp vs. C,e, with E"^ = p"^ + Qq, 
Qo = 270 MeV at various cms energies with predictions of the MLLA iterative 
solution (H) (Qo = 270 MeV, A = 0.01, Kh = 0.45). Data and predictions 
for particle energies E < 1 GeV are shown. 

Fig. 9: The DLA asymptotic prediction ( ^6|) for the scaling function E{(), 
the prediction of the Gaussian Approximation (^Sp and the corresponding 
predictions with fixed in comparison. 

Fig. 10: Test of ^-scaling for the momentum spectra at 14, 58 and 130 
GeV |T0| , |26|Jl^ with average charged multiplicities taken from [^,^,0 . 

Fig. 11: Test of ^-scaling for the energy spectra. Data as in Fig. 10, but the 
particle energy is calculated using the mass Qo = 270 MeV, in comparison 
with theoretical predictions from the limiting spectrum (Qo = 270 MeV) nor- 
malized to the predicted average multiplicity according to eqs. (^) and 
The DLA prediction ( ^6]) at asymptotic cms energy is also shown. 



Fig. 12: Maximum of the rescaled inclusive momentum distribution (* = 
^* /Y as a function of F = In comparison between experimental data at 
various cms energies [ p!0| , p!l| , p!^ , pr^ , p!^j2^j25|j26| , |27|j28| , |29| and theoretical predic- 
tion in MLLA, numerically extracted from the shape of the limiting spectrum 
(solid line) for the cut-off parameter Qo = A = 270 MeV Crosses mark 



the predictions at the cms energies 200 GeV and 500 GeV. The asymptotical 
DLA result C* = 5 is also shown (dashed line). 
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Table 1 
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MARK II 
5.2 



MARK I0 

7.4 
TASSOp 

14. 



TASSOp 
22. 



TASSOpi 
35. 



TASSOp 
44. 



TOPAZ p I 
58. 



3.01±0.3 
(2.10,input) 
3.66±0.37 
(2.73,3.61) 

± 

(3.36,4.21) 
4.63±0.46 
(4.11,4.92) 
7.64±0.59 
(6.73,7.40) 
9.65±0.68 
(8.87,9.44) 
12.21±0.86 

(11.51,11.95) 
13.38±1.05 

(12.96,13.33) 
14.54±0.43 
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1.02±0.02 

(1.14) 
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(1.33) 
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(1.66) 
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(2.12) 
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0.14±0.01 

(0.13) 
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-0.56±0.08 
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(-0.38) 
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(-0.52) 
-0.29±0.18 

(-0.52) 
-0.42±0.17 

(-0.52) 
-0.44±0.18 

(-0.52) 
-0.51±0.09 

(-0.52) 
-0.59±0.09 

(-0.52) 
-0.49±0.15 

(-0.52) 



ALEPH|7| 
DELPHlH] 



L3[B|] 



OPALlll 



LEP-1 (avg) 
DELPHI-7 111 
91.2 



18.81±1.05 
19.17±1.00 
18.74±1.09 
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18.93±0.52 
19.20±0.26 
(input, input) 
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3.29±0.01 
3.33±0.04 
(3.27) 



0.99±0.05 
1.03±0.01 
0.99±0.06 
0.99±0.01 
1.01±0.02 
0.99±0.03 
(1.00) 



-0.39±0.10 
-0.40±0.02 
-0.35±0.13 
-0.36±0.03 
-0.39±0.02 
-0.46±0.05 
(-0.35) 



-0.59±0.32 
-0.59±0.07 
-0.65±0.40 
-0.59±0.09 
-0.59±0.05 
-0.44±0.18 
(-0.52) 



ALEPH 112 



DELPHI lH] 

OPAL 
LEP-1. 5 (avg; 
133 



22.04±0.47 
22.27±0.58 
21.50±0.57 
21.95±0.31 
(22.6,22.5) 



3.52±0.06 
3.47±0.05 
3.51±0.07 
3.495±0.034 
(3.50) 



1.19±0.04 
1.13±0.05 
1.19±0.04 
1.185±0.025 
(1.14) 



-0.37±0.07 
-0.40±0.08 
-0.35±0.07 
-0.365±0.042 
(-0.34) 



-0.62±0.26 
-0.49±0.29 
-0.63±0.25 
-0.59±0.15 
(-0.52) 
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